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Abstract. We develop a theory for filamentary molecular clouds in- 
cluding the effects of ordered magnetic fields, and external pressure. We 
first derive a new virial equation appropriate for filamentary clouds. By 
comparing with observational results collected from the literature, we 
find that the fields are likely helical. Secondly, we construct numerical, 
MHD models of filamentary clouds that agree with the observational con- 
straints. We find that our models produce more realistic density profiles 
r ~ r ^ 18 to ~ 2 than previous models, where the density falls off as r~ 4 . 



1. Introduction 

Most molecular clouds are filamentary structures that are supported by non- 
thermal, MHD turbulence, as well as large scale ordered magnetic fields (eg. 
review, McKee et. al 1993). Nevertheless, virtually all theoretical models assume 
spheroidal geometry. While spheroidal models obviously apply to molecular 
cloud cores, they cannot adequately describe molecular clouds on larger scales. 

Observations suggest that some filamentary clouds may be wrapped by 
helical fields (Heiles 1987; Bally 1987). A few authors have previously modeled 
filamentary clouds with helical fields (cf. Nakamura, Hanawa, & Nakano 1993; 
Hanawa et al., 1993), but the fields in these models simply rescale the Ostriker 
(1964) solution for unmagnetized filaments. They are also unconstrained by 
observational data. 

The role of the external pressure has been almost completely ignored by 
previous models of filamentary clouds. As we discuss in Section |2l real molecular 
clouds are truncated at finite radius by the pressure of the external medium. By 
ignoring the external pressure, most existing models cannot adequately describe 
real molecular clouds. 



2. Surface Pressures on Molecular Clouds 

Molecular clouds and the surrounding atomic gas are dominated by non-thermal 
motions, which result in total pressures that greatly exceed the thermal pressure. 
The total pressure of the atomic gas has been evaluated at the Galactic midplane 
by Boulares and Cox (1990), who find pressures on the order of 10 4 K cm~ 3 . 
However, it is likely that the molecular clouds are surrounded by atomic gas at 
significantly higher pressures. Some are associated with HI clouds at a pressure 
of ~ 10 5 K cm~ 3 (Chromey, Elmegreen, and Elmegreen 1989). Therefore, we 
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conservatively adopt surface pressures in the range of 10 4 5 K cm 3 for our 
analysis. 

3. A New Form of the Virial Equation for Filamentary Molecular 
Clouds 

We consider the virial equilibrium of a long filamentary cloud with an external 
pressure Ps that truncates the cloud at some cylindrical radius Rs, and a general 
magnetic field of helical geometry. By considering only radial equilibrium, we 
derive a new virial equation from the tensor virial theorem: 

where 
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In this equation, Pg and (P) are the surface and average internal pressures, m 
is the mass per unit lengh, m v i r is the "virial" mass per unit length. W and M 
are respectively the gravitational and total magnetic energies per unit length. 
We note that W is independent of the filament radius. 

The total magnetic energy M may be either positive or negative depending 
on whether the poloidal field or the toroidal field dominates the overall energetics 
of the cloud. Thus, the poloidal field component helps to support the cloud 
radially against self-gravity, while the toroidal component works with gravity 
and the external pressure to squeeze the cloud. 

The virial quantities m/m v i r and Pg/{P), in equation [l], may be determined 
from observations; thus, we may easily infer the magnetic parameter A4/|W|. 
We have calculated m/m v i r and Ps/(P) for several filaments, based on obser- 
vational results gathered from the literature. The reader is referred to Fiege & 
Pudritz (1999a) for the full data tables and references. We refer to Figure 1 
of Pudritz and Fiege (1999), in this volume, where we plot contours of A^/|W| 
over these virial parameters. We find that most of the filaments in our sample 
fall in the range 

0.11 < m/m vir < 0.43 
0.012 < P S /(P) < 0.75, (3) 



and that _M/|W| < 0, which is consistent with a helical magnetic field. 
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Figure 1. The density and magnetic structure of our models. The 
dashed lines represent r -4 and r~ 2 profiles for comparison. 

4. Numerical Magnet ost at ic Models 

Our theoretical models involve three parameters, two to describe the mass load- 
ing of the poloidal and toroidal field lines, and a third that specifies the radial 
concentration of the filament. We fully sample our parameter space by a Monte 
Carlo method, in which models are randomly generated and tested for agreement 
with the observational constraints given in equation |3[ In addition to equation 
||, we demand that the magnetic and kinetic energies are nearly in equipartition; 
0.2 < ^? < 5, where M and K are the average magnetic and kinetic energy den- 
sities in the cloud. This has been observationally determined for many clouds 
(Myers & Goodman 1988a,b). 

Figure [l] shows the structure of a few models that are consistent with our 
constraints. We find that the density falls off as r -1 ' 8 to r -2 in the outer regions, 
which is in excellent agreement with recent results by Alves et al. (1998) and 
Lada et al. (1998) for the L977 and IC5146 filamentary molecular clouds. We 
find that the peak poloidal field in our models is always much stronger than the 
peak toroidal field; only the outer regions are dominated by the toroidal field 
component. In this sense, our helical fields are actually quite weakly wrapped. 
These results are robust; a Monte Carlo sampling of our parameter space shows 
that nearly all allowed models have these characteristics. 

5. Stability of Filamentary Molecular Clouds 

We consider the stability of filamentary molecular clouds against axisymmetric 
fragmentation into cores. By solving the linearized equations of MHD and self- 
gravity, we determine how pressure truncation and the individual magnetic field 
components affect the stability of our models. We find that pressure truncation 
has a stabilizing effect, which results from the decreased mass per unit length of 
the filament. Both field components also stabilize clouds against gravitational 
fragmentation. However, we find that sufficiently strong toroidal fields may trig- 
ger MHD-driven instabilities. We refer the reader to Fiege & Pudritz (1999b), 
where we fully explore the stability of our models. 
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6. Discussion 



We have found indirect evidence, based on a very general virial analysis, suggest- 
ing that most filamentary molecular clouds may be wrapped by helical magnetic 
fields. We have also constructed numerical MHD models that agree with the 
observations. These models differ from previous models in that they have much 
shallower density profiles, which are in good agreement with the available ob- 
servational results; we find that the density falls off as r~ 1,8 to r~ 2 , compared 
with the r -4 behaviour of previous models. This behaviour is entirely due to 
the toroidal character of the magnetic field in the outer regions. We also find 
that models with purely poloidal fields have steep density gradients that are not 
allowed by the observations. 
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